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Abstract 

Background: Cervical auscultation (CA) is an affordable, non-invasive technique used 
to observe sounds occurring during swallowing. CA involves swallowing 
characterization via stethoscopes or microphones, while accelerometers can detect 
other vibratory signals. While the effects of fluid viscosity on swallowing accelerometry 
signals is well understood, there are still open questions about these effects on 
swallowing sounds. Therefore, this study investigated the influence of fluids with 
increasing thickness on swallowing sound characteristics. 

Method: We collected swallowing sounds and swallowing accelerometry signals from 
56 healthy participants. Each participant completed five water swallows, five swallows 
of nectar-thick apple juice, and five swallows of honey-thick apple juice. These 
swallows were completed in neutral head and chin-tuck head positions. After 
pre-processing of collected signals, a number of features in time, frequency and 
time-frequency domains were extracted. 

Results: Our numerical analysis demonstrated that significant influence of viscosity 
was found in most of the features. In general, features extracted from swallows in the 
neutral head position were affected more than swallows from the chin-tuck position. 
Furthermore, most of the differences were found between water and fluids with higher 
viscosity. Almost no significant differences were found between swallows involving 
nectar-thick and honey-thick apple juices. Our results also showed that thicker fluids 
had higher acoustic regularity and predictability as demonstrated by the 
information-theoretic features, and a lower frequency content as demonstrated by 
features in the frequency domain. 

Conclusions: According to these results, we can conclude that viscosity of fluids 
should be considered in future investigations involving swallowing sounds. 

Keywords: Swallowing, Swallowing sounds, Viscosity, Signal characteristics 



Introduction 

Dysphagia is a swallowing disorder [1] typically occurring in patients who suffer from 
a variety of neurological conditions (stroke [2], cerebral palsy [3], Parkinsons and other 
neurodegenerative diseases [4]), head and neck cancer and its treatment [5], iatrogenic 
conditions or trauma [6] and other diseases. Dysphagia can also occur due to genetic pre- 
dispositions or congenital craniofacial syndromes [7]. Among the signs and symptoms 
of dysphagia include the subjective sensation of difficulty swallowing food or liquids, 
choking or coughing before, during or after swallowing, or other symptoms caused by 
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impaired clearance of swallowed material into the digestive system, which can cause 
malnutrition [8], dehydration [9], failure of the immune system [10], psycho-social degra- 
dation [11,12] and in general, a decreased quality of life [13]. A major consequence of 
dysphagia is aspiration of food and liquids into the airway past the vocal folds and into 
the respiratory system which leads to airway obstruction, pneumonia, and increased risk 
of mortality resulting from both [14,15]. 

There are several techniques for diagnosing dysphagia using imaging instrumentation. 
The videofluoroscopic swallowing study (VFSS) and the fiber-optic endoscopic evalua- 
tion of swallowing are the currently accepted imaging gold standards [1,16]. These diag- 
nostic methods are available in acute care hospitals, and some rehabilitation centers and 
outpatient clinics, but in some settings such as nursing homes or skilled nursing facilities, 
they are not always readily available and patients need to be scheduled to receive these 
tests at a later date in the acute care hospitals. Furthermore, in some settings, immediate 
performance of imaging studies necessary to definitively diagnose structural or physio- 
logical swallowing disorders leading to increased risk of adverse medical events is not 
possible at a moment s notice, leaving clinicians to use screening methods in an attempt 
to predict likely impairments and manage them while awaiting imaging assessments, 
despite the low precision of screening tests in the identification of impairments. In the 
screening of stroke patients at the acute care setting, a widely accepted standard practice 
in the US, immediate swallowing screening is performed upon immediately admission 
before the patient has had an opportunity to eat or drink or take oral medications, to 
identify likelihood of aspiration because aspiration and its adverse outcomes significantly 
increases morbidity after stroke [17]. Though imaging procedures carry some degree of 
invasiveness such as exposure to radiation and intubation by a fiber-optic endoscope, they 
remain necessary for accurate identification of impairments of swallowing function and 
determination of treatment options that might alleviate impairments or lower the risks 
caused by the impairments. They must also be performed by a trained diagnostic spe- 
cialist. Therefore interest in less invasive screening methods has gained momentum over 
the past 20 years. A non-invasive method of screening for dysphagia known as cervical 
auscultation (CA) has been explored in recent years [18], although its ability to identify 
or predict specific features of dysphagia or guide intervention to alleviate risks asso- 
ciated with dysphagia has not been established [19]. CA usually involves investigating 
signals acquired via stethoscopes or microphones [20,21]. As in all noninvasive screen- 
ing methods, one attraction of CA is its mobility for day-to-day monitoring and a low 
price [21] though its predictive value for identifying important diagnostic signs has yet to 
be established. CA as a tool for screening for dysphagia is still under investigation (e.g., 
[18,21]). 

Previous studies indicated that thicker liquids can reduce the amount of material that 
is aspirated when individuals aspirate thin liquids while swallowing [22] or subjectively 
improve swallowing symptoms in some individuals who have dysphagia with ordinary 
liquids so it would be informative to determine whether the effects of increased fluid vis- 
cosity on swallowing signal characteristics produces useful information that might add 
value to auscultation as a screening method [23,24]. Though there is understanding of 
the effects of increased viscosity on swallowing accelerometry signals (e.g., [25]), the 
effects on swallowing acoustics are more challenging to understand. One challenge is that 
previous studies used microphones of a varying quality to acquire swallowing sounds. 
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In [26], the authors used Sony ECM-C115 microphone with a frequency response from 
50 Hz to 15 kHz to show that duration of the swallow signals are longer for thicker 
fluids. A similar trend was observed by Reynolds et al [27] using an electret micro- 
phone Optimus (Radio-Shack/Tandy Corp, Model 333013), with a nonlinear frequency 
response form 70 Hz to 16 kHz. Other challenges to the usefulness of auscultation in 
dysphagia screening stems from the previously adopted microphones, which were not 
able to capture low frequency components of swallowing sounds. In our recent study 
Dudik JM, Jestrovic I, Luan B, Coyle J, Sejdic E: A comparative analysis of swallowing 
accelerometry and sounds during saliva swallows. [Submitted], we showed that the swal- 
lowing sounds are centered at lower frequencies below 50 Hz and their bandwidth extends 
up to few hundred Hertz. These open challenges prompted us to conduct the current 
investigation. 

In this paper, we sought to investigate the effects of fluids with increased viscosity 
on swallowing sound characteristics. In particular, we examine the signal characteristics 
in time, frequency and time-frequency domains, while participants completed swallows 
in neutral head-neck posture and the head-neck flexion (chin-tuck) position which has 
also been used to manage aspiration in patients with specific biomechanical swallow- 
ing impairments [28-31]. To compare our results with the previous study [25], we also 
simultaneously collected dual-axis swallowing accelerometry signals. 

Methodology 

Data acquisition from participants 

In this study, simultaneous accelerometry and acoustic data were collected from 56 
healthy adults aged 18 to 65 years. All subjects had no previous history of neuro- 
logical diseases, swallowing difficulties and/or cancer of the mouth, neck or brain. 
The study protocol was approved by Institutional Review Board at the University of 
Pittsburgh. 

After signing a consent form and recording information about subjects' height and 
weight, the dual-axis accelerometer (ADXL322, Analog Devices, Norwood, MA, USA) 
and contact microphone (AKG C411L, AKG Acoustics GmbH, Vienna, Austria) were 
attached to the subjects neck with double sided tape. The accelerometer was placed below 
the thyroid cartilage and the microphone was placed below the accelerometer, far enough 
to avoid contact between two sensors, as shown in Figure 1. The experimental procedure 
was divided in to two parts conducted in the same order for all participants. First, partic- 
ipants completed bolus swallows in a neutral-head position, followed by the completion 
of swallows in a chin-tuck position. In both parts, the subject was asked to take five indi- 
vidual swallows of different fluids: water, nectar consistency and honey consistency apple 
juices. Thickened apple juices are commercially available products (Nestle Health Care, 
Inc. Florham Park, NJ, USA). Nectar consistency and honey consistency apple juices are 
classified by the Australian Standard for Texture Modified Foods and Fluids, as Mildly 
Thick-Level 150 for nectar and Moderately Thick-Level 400 for honey- thick. All fluids 
were served chilled (3-5°C) in cups as approximately one bolus per cup. Participants were 
asked to complete the individual swallows of a single bolus at a comfortable pace while 
consuming comfortable volumes of fluids. The volume of bolus was not controlled as 
there are sex based differences in a comfortable bolus size [32]. We intend to investigate 
the effects of specific bolus volumes in future research. 
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Figure 1 Position of accelerometer and microphone. 



The accelerometer was powered with a 3V output power supply (1504 DC/ AC Power 
Supply, B&K Precision Corporation, Yorba Linda, CA, USA). The two axes of the 
accelerometer were positioned in the anterior-posterior (A-P) and superior-inferior (S-I) 
directions. Signals from both axes of the accelerometer were passed through an amplifier 
(P55, Grass Technologies, Warwick, RI, USA), which provided 10 times amplification and 
then were band-pass filtered from 0.1 to 3000 Hz. A microphone was powered by a power 
supply (model B29L, AKG, Vienna, Austria). Swallowing accelerometry signals and swal- 
lowing sounds were sampled at 40 kHz by the LabView program Signal Express (National 
Instruments, Austin, TX, USA) running on a personal computer. All data were saved on 
an external hard drive. 

Pre-processing steps 

First, all collected signals were pre-processed according to previously proposed algo- 
rithms (e.g., [25]). The accelerometer signals were downsampled to 10 kHz in order to 
implement previously proposed approaches [25]. 

All acquired signals were initially filtered with a finite impulse response (FIR) filter to 
annul the effects of the data acquisition equipment. The filters for swallowing accelerom- 
etry signals and swallowing sounds were designed according to the procedure outlined in 
[33] using 18 table-top recordings in a quiet room. 

Next, we removed very low frequency components from the dual-axis accelerome- 
try signals associated with head movements [34]. Since the microphone signal was not 
affected by any head movements, there was no need to perform such an operation for 
these signals. 

Consequently, all signals were denoised with 10-level discrete wavelet decomposi- 
tion using the discrete Meyer wavelet with soft-thresholding using the global denoising 
threshold, T^ en defined as: 
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_ med(\di\)y/2 logn 
den ~ 0.6745 ' W 
where d\ represents wavelet coefficients at the first level, n is length of the signal and med 
is median operator [35]. 

The last pre-processing step was the segmentation of signals carried out according to 
the sequential fuzzy c-means algorithm designed for dual-axis accelerometry signals [36] . 
All segmentation results were verified visually, if any of them were incorrect, swallows 
were segmented manually. Swallows which could not be segmented were excluded from 
the study (less than 5% were excluded). The time instances identified in this process as 
the beginning and the end of each swallow were then used to segment the microphone 
signals. 

After the completion of all pre-processing steps, the features outlined in the next 
subsection were extracted from each swallow. 



Feature extraction 

Each swallowing sound could be represented as a discrete time series, M = {mi, 
m2t . . .,m n }. Different signal features can be used to describe swallowing characteristics, 
and we summarize below the features considered in this study. The same set of fea- 
tures was considered for both swallowing sounds and dual-axis swallowing accelerometry 
signals. 

Time domain features 

• The mean (average) value of a signal represents unbiased estimation of the amplitude 
of the signal. An equation for calculating the mean value is given as 

1 n 

Mm = - m. (2) 

i=i 

• The standard deviation is a measure of variation from the mean value. It can be 
obtained as 



a = 

N 



i n 

V (mi - fi m ) 2 . 

n — l f—* 



(3) 



The skewness represents symmetry of a distribution of the signal [25] . It can be 
calculated as, 



(4) 



The kurtosis is a measure of the "peakedness" of the probability distribution of a 
variable. For a high value of kurtosis, the distribution is sharp and narrow, with heavy 
tails. A low kurtosis value indicated a flat distribution peak and thin tails. Kurtosis is 
calculated as 

The entropy rate [37,38] quantifies the extent of regularity in a signal. It provides 
important information about swallows as an random process. Entropy rate is 
calculated in several steps. First, a signal M should be normalized to zero mean and 
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unit variance. The normalized M is then quantized to 10 equally spaced levels. Those 
10 levels are ranged from minimum to maximum and marked with integer numbers 
from 0 to 9. Then the quantized signal M = {mi, m<i> . . .,rh n }, with U consecutive 
points is coded as 

Si = m+u-i ' lO^ 7 " 1 + • • • + m • 10°, (6) 

where i = 1, 2, . . . , n — U + 1, and S\ = {si, $2> • • • Sn-u+i\ are coded integers. 
Because of the 10 quantization levels, 10 is used as a base. Using the Shannon entropy 
formula, the entropy is estimated as 

lO^" 1 

E{U) = -J2 p SuO<) ' lnPs u (k)> (7) 

k=l 

where Ps u is probability of observing k in S u , approximated by the corresponding 
sample frequency. The entropy is then normalized using following formula 

- E(U) -E{U- 1) +£(1) -a 

NE(U) = — K - J — 9 (8) 

where a is the percentage of the coded integers in Si that occurred only once. Finally, 
the regulatory index as a measure of the entropy rate is calculated as 

p = 1 - minNEXU). (9) 

p takes value from 0 to 1, where for regulatory index is equal to 1 indicates maximum 
of regularity, while value of 0 represents maximum of randomness. 

• The Lempel-Ziv complexity (L-Z) [39] provides information about predictability of 
the signal. To compute the L-Z complexity, a signal M should be first quantized into 
100 equally spaced levels. Then this 100 levels are ranged from minimum to 
maximum values. In the next step, the quantized signal A\ = {a\, a2, • . • , a n } was 
decomposed in L different blocks of the length / — ; + 1, so that A\ = {\J/i, \J/2, . . . , 
\j/ n }. Blocks are defined as 

vl/ = A\ = {aj, aj+i, . . . , */}, 1 < j < I < n (10) 

The first block is equal to the first element of the quantized signal. Other blocks are 
defined as 

^ +1 =A^ + 1 1 ,m>l,mGZ+ (11) 
where h m is ending index for \j/ m . Finally, the L-Z complexity is calculated as 

LZ = Il0gl °°* (12) 
n 

Frequency domain features 

• The peak frequency of a signal is defined as 

^ = argmax|F M (f)| 2 , (13) 

feWmax] 

where f max is the highest available frequency in a signal and Fm represents the 
Fourier transform of a signal. 
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• The centroid frequency indicates position of the center of mass in the signal in the 
frequency domain [33]. For the signal, M, it is estimated as 

fmax 

f f\F M (f)\ 2 df 

fc = -7 • (14) 

Jmax 

f \F M (f)\ 2 df 
o 

• Bandwidth represents spectral spread and it is defined as 



BW = 



fmax _ 

/ (f -fc) \FM(f)\ 2 df 
_0 

fmax 

f \FM(f)\ 2 df 

N 0 



(15) 



Time frequency domain feature 

• The relative energy was computed using a 10-level discrete wavelet decomposition of 
the signal with the Meyer wavelet [25,40-42]. The energy at each decomposition level 
is computed using the Euclidean norm of decomposition coefficient vectors: 

E ai0 = \\a 10 1 \ 2 > (16) 

E di = \\di || 2 , (17) 
where a\Q is the approximation signal and d{ is detail signal. The total energy was 
calculated as 

10 

i=l 

Finally, percent of relative energy contribution from each decomposition level was 
computed as 

Et ai0 = ^ x 100%, (19) 



E td . = x 100%, (20) 



E T 

for i = 1, 2, . . . , 10. 

• Wavelet entropy describes the information distribution in the time-frequency 
domain. Wavelet entropy was computed using 10-level wavelet decomposition and 
relative energy computed above, with following formula: 

WE = • logfi^ -Y^L. log 2 ^- , (21) 

100 52 100 ^100 52 100 

i—l 

Data analysis 

The statistical differences between different conditions were tested using a non- 
parametric statistical hypothesis test, Wilcoxon rank-sum test [43]. 



Results 

Results of the feature extraction process are presented as a mean value =b standard devi- 
ation. We analyzed 271 water swallows in neutral position and 274 in chin- tuck position, 
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277 nectar-thick apple juice in neutral position and 275 in chin-tuck position, and 273 
honey-thick apple juice swallows in neutral position and 273 in the chin tuck position. 

Time domain features results 

Table 1 summarizes the time domain features from the swallowing sounds. The results 
showed that standard deviation (a), skewness (v) and kurtosis (m) were not significantly 
different between the control condition (water) and the thickened liquid conditions in 
the chin- tuck position (p > 0.05, zval < —0.49, ranksum < 76218). For the swallows 
in the neutral position, pairwise comparison between water and nectar-thick apple juice 
revealed statistically significant differences for standard deviation (p = 0.03, zval = 2.23, 
ranksum = 81398) and skewness (p — 0.01, zval = —2.53, ranksum = 72315). The skew- 
ness was significantly different between water and honey-thick apple juice (p << 0.01, 
zval = —3.73, ranksum = 69545) as well as the kurtosis (p = 0.02, zval = 2.51, 
ranksum = 765150). Next, we observed significantly higher entropy rates (p) for nectar- 
thick and honey-thick fluids in comparison to water for both head positions (p << 0.01, 
zval < —4.09, ranksum < 70304). However, the L-Z complexity had statistically the 
highest values for water swallows for both head maneuvers (p < 0.05, zval < 3.22, 
ranksum < 805970). 

Table 2 summarizes the results for the swallowing accelerometry signals. The results 
showed that in the A-P direction of the accelerometer signal, standard deviation and 
kurtosis in the chin-tuck position were not affected by the fluid viscosity (p > 0.05, 
zval < 2.69, ranksum < 82959). Water swallows in the neutral position had the statis- 
tically highest values for standard deviation (p < 0.01, zval < 3.67, ranksum < 83135) 
and the lowest values for kurtosis (p < 0.03, zval < 3.95, ranksum < 84631). The skew- 
ness was statistically different between nectar-thick and honey-thick apple juice in neutral 
position (p = 0.03, zval = — 2.01, ranksum = 74794), and between water and honey- 
thick apple juice in chin-tuck position (p = 0.03, zval = —2.48, ranksum = 73126). 
Furthermore, water swallows had statistically the lowest values for entropy rate (p < 
0.05, zval < —5.01, ranksum < 68555) and the highest values for the L-Z complex- 
ity (p << 0.01, zval < 6012, ranksum < 834210) in comparison to other two fluids 
in both head positions. Also, a pairwise comparison between nectar-thick and honey- 
thick swallows found significant differences for entropy rate (p = 0.01, zval = —2.03, 
ranksum = 74383) and L-Z complexity (p = 0.03, zval = —2.87 , ranksum = 83728) in 
the head chin-tuck position. 



Table 1 Time domain features for swallowing sounds 







Neutral position 






Chin-tuck position 




Feature 


Water 


Nectar-thick 
apple juice 


Honey-thick 
apple juice 


Water 


Nectar-thick 
apple juice 


Honey-thick 
apple juice 


(7 


0.54 ± 0.03 


0.42 ± 0.02 


0.54 ± 0.03 


0.54 ± 0.02 


0.54 ± 0.02 


0.54 ± 0.02 


V 


-1.34 ±0.22 


-0.80 ± 0.20 


-1.04 ±0.34 


-1.53 ±0.41 


-2.1 9 ±0.59 


-0.69 ± 0.43 


m 


92.5 ±17.1 


96.1 ± 16.7 


173 ±43.1 


157 ±37.5 


300 ±57.7 


227 ±41 .6 


P* 


98.7 ± 0.04 


99.0 ± 0.04 


99.1 ±0.06 


98.1 ±0.14 


98.5 ±0.10 


98.7 ± 0.05 


LZ* 


6.14 =b 0.15 


5.78 ±0.16 


5.61 ±0.18 


7.45 ± 0.29 


6.39 ± 0.26 


5.98 ± 0.20 



* denotes multiplication by 10" 2 . 
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Table 2 Time domain features for swallowing accelerometry signals 







INcUlldl pualliuil 






V.I1 III LUV.lV [JUblllLMI 




Feature 


Water 


Mortar-thirls 

apple juice 


nui icy ii iiv.iv 

apple juice 


Water 


iNcvAdl llllv.lv 

apple juice 


nuiitryiiiiv.lv 

apple juice 


er* A-P 


1.39±0.05 


1 .1 6d=0.03 


0.39±0.02 


1 .39±0.04 


1 .39±0.04 


1 .39±0.04 


er* S-l 


1.1 1 =b0.06 


0.96±0.03 


1.16±0.05 


1.16±0.05 


1.16±0.04 


1.16±0.05 


v A-P 


-0.73±0.22 


-1.39±0.23 


-0.74±0.21 


-2.31 ±0.43 


-2.24±0.49 


-1.31 ±0.42 


v S-l 


0.28±0.32 


0.14±0.37 


-0.49±0.39 


-0.13±0.31 


-0.69±0.29 


-0.54±0.37 


m A-P 


64.5±12.8 


62.7±16.7 


64.1 ±13.6 


173±30.5 


193±42.1 


183±33.6 


m S-l 


81.8=b1 7.0 


121 ±28.1 


118±32.2 


96.9±21.2 


193±21.5 


145±22.6 


p* A-P 


98.8±0.04 


99.1 ±0.02 


99.1 ±0.04 


98.5±0.07 


98.8±0.06 


99.1 ±0.04 


p* S-l 


99.1 ±0.03 


99.2±0.02 


99.2±0.03 


98.5±0.08 


98.8±0.04 


98.9±0.04 


LZ* A-P 


5.46±0.12 


4.97±0.12 


4.92±0.14 


6.26±0.19 


5.44±0.17 


4.83±0.14 


LZ* S-l 


6.36±0.14 


6.21±0.15 


6.31±0.16 


7.17±0.22 


6.42±0.21 


5.91 ±0.18 



* denotes multiplication by 10~ 2 . 



In the S-I direction, the fluid thickness did not have influence on L-Z complexity in 
the head neutral, and standard deviation and kurtosis in the chin-tuck position (p < 
0.05, zval < 2.01, ranksum < 81661). For skewness, nectar swallows showed a signifi- 
cant statistical difference in neutral position (p < 0.02, zval < 2.49, ranksum < 83601), 
while in chin-tuck position water swallows has the lowest value (p < 0.02, zval < 
2.45, ranksum < 83338). The standard deviation was statistically different between water 
and nectar-thick (p — 0.02, zval — 1.41, ranksum — 80110) as well as kurtosis between 
water and honey-thick (p = 0.01, zval = 2.95, ranksum = 82742). Additionally, the 
entropy rate is observed to be significantly lower in water swallows than in the other 
two stimuli in both head position (p << 0.01, zval < —3.42, ranksum < 71586). 
Water swallows showed a significantly higher value for the L-Z complexity in the chin- 
tuck position (p < 0.05, zval < 4.37 , ranksum < 86156), while a pairwise comparison 
between nectar-thick and honey-thick apple juices showed a difference for the entropy 
rate (p = 0.02, zval = -1.71, ranksum = 74989). 

Also, we compared the extracted features between two accelerometer axes. Kurtosis in 
both head positions did not exhibit a significant statistical difference (p > 0.05, zval < 
0.91, ranksum < 78599). The standard deviation in the neutral head position and 
skewness in the chin-tuck position showed statistically significant differences between 
swallows for all stimuli (p << 0.01, zval < 1 71% ranksum < 92505). In the neutral 
position, skewness was significantly different between water and nectar-thick swallows 
(p << 0.01, zval < —4.64, ranksum < 69137), while the standard deviation showed 
a significant difference for nectar-thick swallows in chin-tuck position (p << 0.01, 
zval < —3.07, ranksum < 74432). The L-Z complexity and the entropy rate were also 
significantly different for all stimuli in both head positions (p << 0.01, zval < 1.31, 
ranksum < 715650). 

Frequency domain features results 

Table 3 summarizes the values of the considered frequency features for swallowing 
sounds. The centroid frequency (f c ) and the bandwidth (BW) were not affected by the 
fluid viscosity in the chin-tick position (p > 0.05, zval < 0.62, ranksum < 79457), 
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Table 3 Frequency domain features for swallowing sounds 







Neutral position 






Chin-tuck position 




Feature 


Water 


Nectar-thick 


Honey-thick 


Water 


Nectar-thick 


Honey-thick 






apple juice 


apple juice 




apple juice 


apple juice 


f P 


26.6±4.93 


16.7zb1.96 


8.68=b1 .69 


24.3±3.82 


17.9±2.19 


13.5±1.71 


fc 


446±45.4 


464±51.6 


493±65.7 


739±66.1 


802±69.5 


767±73.4 


BW 


759±46.3 


736±60.3 


725±61.1 


11 61 ±68.5 


1269±72.6 


1236±71.8 



while the peak frequency (f p ) had significantly higher values for water swallows in the 
chin-tuck position (p < 0.04, zval < 4.46, ranksum < 82550). In the neutral head 
position, the peak frequency was significantly higher for water swallows than for honey- 
thick swallows (p = 0.01, zval = 2.29, ranksum < 80912), while simultaneously the 
water swallows had significantly smaller bandwidth values than the honey-thick swallows 
(p = 0.02, zval = 2.49, ranksum = 81282). The water swallows also had the small- 
est values for the centroid frequency in comparison to the other two types of swallows 
(p << 0.01, zval < 3.81, ranksum < 817740). 

The centroid frequency and bandwidth of the swallowing accelerometry signal in 
the A-P direction was not affected by fluid viscosity in the chin-tuck position (p > 
0.05, zval < 0.93, ranksum < 79973). However, in the A-P direction the centroid fre- 
quency and bandwidth has significantly higher value for water swallows in the neutral 
position (p « 0.01, zval < 4.81, ranksum < 86252). In the same direction, a pair- 
wise comparison between water and honey-thick apple juice for the peak frequency 
showed differences in neutral position (p = 0.006, zval = 2.74, ranksum < 823350), 
while in the chin-tuck position honey-thick swallows had statistically the lowest value 
(p < 0.02, zval < 4.33, ranksum < 86078). 

In the S-I direction, fluids did not impose any statistical differences on the cen- 
troid frequency in chin-tuck position, nor or the bandwidth in the neutral position 
(p > 0.05, zval < 1.11, ranksum < 79528). The peak frequency was significantly dif- 
ferent only between water and nectar-thick swallows in the head-neutral position (p = 
0.03, zval = 2.06, ranksum = 81347), and between water and honey-thick swallows in 
the chin-tuck position (p < 0.01, zval = 2.85, ranksum = 83258). However, also in the 
S-I direction, the centroid frequency exhibited significant differences between water and 
honey- thick swallows in the neutral head position (p = 0.01, zval = 2.41, ranksum = 
81726), while water swallows had smaller bandwidth values than the nectar-thick 
swallows in the chin- tuck head position (p = 0.02, zval = —2.2047 , ranksum = 
73910). Table 4 summarizes frequency characteristics for swallowing accelerometry 
signals. 

While comparing statistical differences between the A-P and S-I directions, we found 
significant differences for the peak frequency for all stimuli in both head positions 
(p << 0.01, zval < —9.53, ranksum < 598280). Furthermore, the centroid frequency 
is different for nectar-thick and honey-thick swallows in both head positions (p < 0.03, 
zval < — 2.16, ranksum < 74720). The bandwidth was significantly different between 
the two directions for all stimuli (p < 0.04, zval < 5.29, ranksum < 86842) in the neu- 
tral head position. Lastly, the bandwidth was significantly smaller for the S-I direction for 
water swallows in the chin-tuck position (p = 0.02, zval = 2.33, ranksum < 82434). 
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Table 4 Frequency domain features for swallowing accelerometery signals 
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244±17.9 


243±17.6 


BWS-\ 


94.8±9.89 


89.7±9.23 


85.8±11.1 


174±13.3 
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Time-frequency domain feature 

The relative energy decompositions are presented in Figures 2, 3, 4, while the wavelet 
entropy results for both swallowing sounds and accelerometry signals are summarized in 
Table 5. 

The wavelet analysis of the swallows showed that the viscosity of fluids had a major 
impact on the time-frequency structures of these signals. Let us first consider the swal- 
lowing sounds. From Figure 2, it is obvious that majority of the energy is concentrated on 
the first alO level for both head maneuvers. Levels dlO and d9 in the neutral head position, 
and d8, d7, d6, d5, and dl in the chin-tuck position were not affected by viscosity of the 
fluids (p > 0.05, zval < 1.77, ranksum < 82367). In both head positions, water swallows 
had the statistically lowest value in the alO level. However, water swallows had a higher 
energy concentration than the other two stimuli in the most of higher frequency levels 
(d8, d7, d6, d5, d4, d3,and dl (p < 0.04, zval < 5.18, ranksum < 86369)) in the neutral 
head position. Also, nectar swallows were statistically different from other stimuli for lev- 
els d4, d3 and d2 (p < 0.03, zval < 5.016, ranksum < 86038). In chin-tuck head position, 
nectar swallows are shown to have statistical difference from other fluids in levels alO and 
d3 (p < 0.01, zval < 3.11, ranksum < 84170), while water swallow has the lowest value 
at level dlO (p < 0.01, zval < —2.46, ranksum < 73151). A pairwise comparison between 
water and honey-thick apple juice revealed significant differences for levels d4 and d2 
(p << 0.01, zval < 3.11, ranksum < 83755), while water and nectar-thick apple juice 
were significantly different for the level d9 (p = 0.01, zval = —2.41, ranksum = 73515). 
Lastly, the wavelet entropy (WE) had a smaller value for the fluids with higher viscos- 
ity in the neutral head position (p < 0.02, zval < 5.14, ranksum < 86278), while in the 
chin-tuck position, nectar swallows exhibited a significant difference from the other two 
swallows (p < 0.02, zval < 2.64, ranksum < 83295). 
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Figure 2 Mean relative energy per decomposition band for swallowing sounds. 
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Figure 3 Mean relative energy per decomposition band for swallowing accelerometry signals in the 
A-P direction. 



Contrary to the previous study on swallowing accelerometry signals [44], a significant 
influence of fluid viscosity was noticed on the swallowing accelerometry signals from 
both directions. First, let us consider the relative energy decomposition of the swallow- 
ing accelerometry signals in the A-P direction. Similar to the swallowing sounds, most 
of the energy is concentrated in the alO level for all fluids. Additionally, water swal- 
lows have the statistically lowest energy concentration in the alO level (p << 0.01, 
zval < —1.17, ranksum < 755616), which was not the case at higher frequencies, 
where water swallows had mostly higher energy concentration for both head maneu- 
vers and both axes. The results for the A-P direction showed that the dlO and d9 levels 
in the neutral position and most of the levels in chin tuck position were not affected 
with viscosity of fluids. In the neutral head position, all stimuli showed a significant 
difference in the levels alO, d4 and d3, while water swallows exhibited higher energy con- 
centrations in the d8, d7, d6, and d5 levels (p < 0.01, zval < 5062, ranksum < 88134). 
Nectar-thick apple juice swallows revealed a significant difference in the d2 level 
(p < 0.03, zval = 0.41, ranksum = 78187) for the neutral head position. In the chin-tuck 
head position, water swallows showed significant difference in level dl (p < 0.05, zval < 
5.62, ranksum < 88134), while a pairwise comparison between water and nectar-thick 
showed a significant difference in the dlO level (p = 0.01, zval = 1.69, ranksum < 74189). 
Lastly, in the A-P direction, the wavelet entropy had a significantly lower value for flu- 
ids with higher viscosity in the neutral position (p < 0.02, zval < 5.14, ranksum < 
86278). The wavelet entropy was not affected by viscosity in the chin-tuck position 
(p > 0.05, zval < 2.64, ranksum < 83295). 
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Figure 4 Mean relative energy per decomposition band for swallowing accelerometry signals in the 
S-l direction. 
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Table 5 Wavelet entropies for swallowing sounds and accelerometry signals 
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Honey-thick 


Water 


Nectar-thick 


Honey-thick 






apple juice 


apple juice 




apple juice 


apple juice 


WE 


1.81 ±0.04 


1.65±0.04 


1.51 ±0.04 


1.67±0.04 


1.69±0.05 


1.51 ±0.05 


WE A-P 


1.78±0.04 


1.55±0.04 


1.39±0.03 


1.71 ±0.04 


1.65±0.04 


1.65±0.04 


WE S-l 


1.91 ±0.03 


1.81 ±0.03 


1.79±0.03 


1.87±0.04 


1.91 ±0.04 


1 .96±0.04 



In the S-I direction, levels d9, d6, and d2 in the neutral head position, and most of the 
levels in chin-tuck position did not show a significant statistical difference between stim- 
uli. Water swallows were significantly different from other fluids in the alO, dlO, d8, d7, 
d4 and d3 levels in the neutral position (p < 0.04, zval < 3.03, ranksum < 84050) , and 
in the d9 level in the chin-tuck position (p << 0.01, zval < —2.53, ranksum < 73018). A 
pairwise comparison between water and honey-thick apple juice exhibited significant dif- 
ferences for the level d5 and dl (p < 0.01, zval < 2.66, ranksum < 82183) in the neutral 
position and for the level dlO (p = 0.01, zval = —2.54, ranksum = 72315) in the chin- 
tuck position. A pairwise comparison between water and honey- thick apple juice showed 
a significant difference in level dlO (p = 0.01, zval = 2.55, ranksum = 72995) in chin- 
tuck head position, while pairwise between nectar-thick and honey-thick apple juice in 
level dl (p < 0.01, zval = — 2.61, ranksum = 73823) in neutral position showed differ- 
ence. Also, the wavelet entropy had statistically the highest value for water swallows in 
the S-I direction (p < 0.02, zval < 2.73, ranksum < 82639). 

The relative energy distribution between the two axes were significantly different 
between each other. Levels alO, dlO, d9, d8, d7, d5, d3 and d2 in the neutral position 
and levels alO, d4, d3 and d2 in the chin-tuck position showed difference between two 
axes for all three stimuli (p << 0.01, zval < 8.91, ranksum < 94965). Furthermore, 
swallows based on nectar-thick and honey-thick apple juices were also different between 
axes for the dl level in the neutral head position and for the d5 level in the chin-tuck 
position (p < 0.01, zval < —2.55, ranksum < 72576). The relative energy distribution 
for water swallows was significantly different between two axes when considering the d5 
level in the neutral position, and the levels d9, d8, d7 and dl in the chin-tuck position 
(p < 0.01, zval < 3.74, ranksum < 85110). However, the d4 level in the neutral head posi- 
tion and the levels dlO and d6 in the chin-tuck position were not significantly different 
between two axes (p > 0.05, zval < 0.33, ranksum < 78820). 

Discussion 

Time domain features 

Our results suggest that the time domain features for swallowing sounds are not different 
between nectar-thick and honey-thick fluids, while the water swallows had significantly 
different features from the other two fluids. These results imply that the difference in 
viscosity between nectar and honey have a limited effect on the extracted time domain 
features. 

For the swallowing sounds, the negative value for skewness indicates that the probability 
distribution of amplitudes are mostly concentrated on the right side (i.e., stronger/louder 
amplitude values). Larger negative skewness values for swallows in chin- tuck position 
denote that swallows have larger (louder) amplitude values in the chin-tuck position than 
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in the neutral position but viscosity did not affect amplitude. Also, kurtosis tends to 
be higher for higher viscosity fluids. Since kurtosis is a measure of "peakedness" of the 
amplitude probability distribution, the results imply that lower viscosity swallows would 
contain more variant amplitudes in the sound signal [26] . Clinically this result indicates 
that detection of varying viscosities of swallowed fluid might be possible with auscultation 
[23]. 

The entropy rate and the L-Z complexity for swallowing sounds were also influenced by 
viscosity of the fluid. According to Table 1, the mean value for the entropy rate is higher 
when viscosity increases, which implies that regularity of the signal is higher for more vis- 
cous fluids [37,38]. Similarly, a higher value for the L-Z complexity means that swallowing 
sounds are more complex and more unpredictable [45,46]. This is in agreement with pre- 
vious studies of CA that have indicated large amounts of signal variability from subject to 
subject and swallow to swallow. From the Table 1, it is obvious that more viscous fluids 
have a lower mean value of the L-Z complexity, which implies that the signal complexity is 
lower for such fluids. The same results were provided by a previous study of the influence 
of viscosity on the accelerometer signal [44] where is implied that higher viscosity fluids 
tends to behave by better defined patterns. These findings indicate that further research 
into the specific characteristics of swallow sounds under various viscosity, posture, and 
other conditions, needs to be elucidated before auscultation will have more clinical value. 

Swallowing accelerometry signals followed similar trends for the entropy rate and the 
L-Z complexity as shown in Table 2. These results confirm the findings from the previ- 
ous study [44], which showed that regularity and predictability is higher for more viscous 
fluids. Also in the previous study, nectar- thick and honey- thick swallows had smaller 
negative values for skewness in the A-P direction. We confirmed the previous results for 
swallows in the chin-tuck position, but failed to confirm this trend for swallows in the 
neutral position. 

Frequency domain features 

As shown in Table 3, thicker fluids yielded swallowing sounds with lower peak frequen- 
cies, which is already proven by a previous study about acoustic nature of normal swallows 
[47]. A similar trend was observed for swallowing accelerometry signals as well. 

Comparing values for swallowing accelerometry signals from Table 4 with those val- 
ues for swallowing sounds from Table 3, it can be concluded that swallowing sounds have 
much higher frequency content than the swallowing accelerometry signals. However, we 
observed similar trends for features extracted from these two types of signals. Bandwidth 
tends to be lower for higher viscosity fluids, which suggests that the more viscous fluids 
required more time for completion of the swallow [23] . The mean value of the centroid 
frequency for swallowing sounds is not dependent on viscosity, which implies that viscos- 
ity does not affect significantly spectral measure [48], which has also been observed for 
the accelerometer signal. 

Time-frequency domain features 

The time-frequency decomposition of swallowing sounds showed that most of the signal 
energy is concentrated at lower frequencies, as was expected based on the frequency anal- 
ysis of swallowing sounds. Thicker fluids have more energy on the first, lowest frequency 
level, since higher viscosity liquids produce a lower swallowing frequency [23]. Clinically, 
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our results can be attributed to the increased total swallow duration, as previous studies 
showed that the oral and pharyngeal swallow durations increase when subjects swallow 
higher viscosity fluids [49,50]. We consider the wavelet entropy to describe spread of 
energy. According to Table 5, the mean value of wavelet entropy tends to be lower for 
higher viscosity fluids, because the energy concentration is higher at the first level for 
thicker fluids. 

Similarly to the swallowing sounds, most of the energy from the accelerometry signals is 
concentrated at lowest frequency level (alO) for all stimuli. Also, the mean value of relative 
energy in alO level tends to be higer for thicker fluids. These findings explain results for 
mean value of wavelet entropy which tends to be lower for thicker fluids. 

Remarks 

According to results, more differences are observed for features in the neutral than in the 
chin-tuck head position for both swallowing sound and swallowing accelerometry sig- 
nals. Furthermore, this study showed more statistical difference for a greater number of 
features extracted from swallowing accelerometry signals than the previous study [44]. 
In the previous study, most of the statistical differences were based on time domain fea- 
tures [44]. It should be mentioned that the previous study only considered data from 17 
participants. 

Limitations and strengths of the present study 

In this study, swallowing conditions have been administered to the subjects in a specific 
order (water, nectar-thick, honey-thick) implying that we cannot rule out the possibility 
that the order of presentation influenced the results. Also, no inference regarding swal- 
lowing physiology can be made from the results of this study as simultaneous imaging 
was not performed. Future research in this area could compensate for these limitations 
by including simultaneously acquired images and randomizing the order of presentation. 
However, this study has contributed to the general knowledge regarding the usefulness 
of CA as a screening method, as we need to clearly understand if there is any more value 
to CA than was previously reported. Future research should also focus on combining CA 
and swallowing accelerometry in a concurrent design (with imaging). The goal would be 
to determine if the detection accuracy of swallowing physiological impairments increases 
by combining these two sensors, or a higher accuracy is achieved by considering sensors 
independently. Also, such studies would enables us to understand the detection accuracy 
of specific physiologic events of these sensors compared to other screening methods. 

Conclusions 

In this paper, we analyzed the effects of fluid viscosity on swallowing sounds in the normal 
and chin-tuck head positions. Swallowing sounds were collected from 56 healthy partici- 
pants, and signal features were extracted from these sounds. Our analysis yielded several 
important conclusions. First, swallowing sounds contained lower frequency components 
than previously reported. Second, fluid viscosity greatly influenced some of the extracted 
features, especially in the frequency and time-frequency domains. Third, most of the time 
domain features exhibited differences between water and fluids with higher viscosity (i.e., 
nectar- thick and honey- thick fluids). The time domain differences were not dominant 
between nectar-thick and honey-thick fluids. 



Jestrovic et al. BioMedical Engineering Online 201 3, 1 2:90 
http://www.biomedical-engineering-online.eom/content/1 2/1 /90 



Page16of 17 



Abbreviations 

CA: Cervical auscultation; VFSS: Videofluoroscopic swallowing study; FIR: Finite impulse response; a: Standard deviation; 
v: Skewness; w. Kurtosis; p: Entropy rate; L-Z: Lempel-Ziv complexity; f c : Centroid frequency; f p : Peak frequency; BW: 
Bandwidth; WE: Wavelet entropy. 

Competing interests 

The authors declare that they have no competing interests. 
Authors' contributions 

IJ, JMD and BL carried out the data collection. IJ conducted all mathematical computations and statistical tests, and 
wrote the entire manuscript. JLC provided clinical insights regarding swallowing and critically revised the manuscript. ES 
designed the data collection protocol, suggested the main objective of the manuscript, provided supervision, and 
critically revised the manuscript. All authors read and approved the final manuscript. 

Acknowledgements 

The authors would like to thank the participants for their time. 
Author details 

1 Department of Electrical and Computer Engineering, Swanson School of Engineering, University of Pittsburgh, 
Pittsburgh, PA, USA. department of Communication Science and Disorders, School of Health and Rehabilitation 
Sciences, University of Pittsburgh, Pittsburgh, PA, USA. 

Received: 22 May 201 3 Accepted: 30 August 201 3 
Published: 1 0 September 201 3 

References 

1 . Logemann JA: Evaluation and treatment of swallowing disorders. Austin, TX: PRO-ED Inc; 1 998. 

2. Gottlieb D, Kipnis M, Sister E, Vardi Y, Brill S: Validation of the 50 ml3 drinking test for evaluation of post-stroke 
dysphagia. Disabil Rehabil 1996, 18(1 0):529-532. 

3. Rogers B, Arvedson J, Buck G, Smart P, Msall M: Characteristics of dysphagia in children with cerebral palsy. 

Dysphagia 1 994, 9:69-73. 

4. Murray J: Manual of dysphagia assessment in adults. Singular Publishing Group: Sad Diego, CA, USA; 1999. 

5. Gillespie MB, Brodsky MB, Day TA, Sharma AK, Lee FS, Martin-Harris B: Laryngeal penetration and aspiration 
during swallowing after the treatment of advanced oropharyngeal cancer. Arch Otolaryngol - Head Neck Surg 
2005,1 31 (7):61 5-61 9. 

6. Halper AS, Cherney LR, Cichowski K, Zhang M: Dysphagia after head trauma: the effect of 
cognitive-communicative impairments on functional outcomes. J Head Trauma Rehabil 1 999, 1 4(5):486-496. 

7. Wilson GN, Oliver WJ: Further delineation of the G syndrome: a manageable genetic cause of infantile 
dysphagia. J Med Genet 1 988, 25(3):1 57-1 63. 

8. Curran J: Nutritional considerations. In Dysphagia: diagnosis and management. Edited by Groher M. Boston: 
Butterworth-Heinemann; 1 992:255-266. 

9. Smithard D, O'neill P, Park C, Morris J, Wyatt R, England R, Martin D: Complications and outcome after acute 
stroke does dysphagia matter? Stroke 1 996, 27(7): 1 200-1 204. 

1 0. Curran J, Groher ME: Development and dissemination of an aspiration risk reduction diet. Dysphagia 1 990, 
5:6-12. 

1 1 . Riensche L, Lang K: Treatment of swallowing disorders through a multidisciplinary team approach. Educ 
Gerontol: /nrQ 1992, 18(3):277-284. 

1 2. Ekberg 0, Hamdy S, Woisard V, Wuttge-Hannig A, Ortega P: Social and psychological burden of dysphagia: its 
impact on diagnosis and treatment. Dysphagia 2002, 1 7(2):1 39-1 46. 

1 3. McHorney CA, Bricker DE, Robbins J, Kramer AE, Rosenbek JC, Chignell KA: The SWAL-QOL outcomes tool for 
oropharyngeal dysphagiain adults: II. item reduction and preliminary scaling. Dysphagia 2000, 1 5(3):1 22-1 33. 

14. Marik PE: Aspiration pneumonitis and aspiration pneumonia. N Engl J Med 2001 , 344(9):665-671 . 

1 5. Ding R, Logemann JA: Pneumonia in stroke patients: a retrospective study. Dysphagia 2000, 1 5(2):51 -57. 

1 6. Tabaee A, Johnson PE, Gartner CJ, Kalwerisky K, Desloge RB, Stewart MG: Patient-controlled comparison of 
flexible endoscopic evaluation of swallowing with sensory testing (FEESST) and videofluoroscopy. 
Laryngoscope 2006, 11 6(5):821 -825. 

1 7. Katzan I, Cebul R, Husak S, Dawson N, Baker D: The effect of pneumonia on mortality among patients 
hospitalized for acute stroke. Neurology 2003, 60(4):620-625. 

1 8. Leslie P, Drinnan M, Finn P, Ford G, Wilson J: Reliability and validity of cervical auscultation: a controlled 
comparison using videofluoroscopy. Dysphagia 2004, 19(4):23 1-240. 

1 9. Leslie P, Drinnan MJ, Zammit-Maempel I, Coyle JL, Ford GA, Wilson JA: Cervical auscultation synchronized with 
images from endoscopy swallow evaluations. Dysphagia 2007, 22(4):290-298. 

20. Cichero J, Murdoch B: Acoustic signature of the normal swallow: characterization by age, gender, and bolus 
volume. Ann Otol Rhinol Laryngol 2002, 1 1 1 (7 Pt 1 ):623-635. 

21 . Zenner P, Losinski D, Mills R: Using cervical auscultation in the clinical dysphagia examination in long-term 
care. Dysphagia 1 995, 1 0:27-31 . 

22. Logemann JA, Gensler G, Robbins J, Lindblad AS, Brandt D, Hind JA, Kosek S, Dikeman K, Kazandjian M, Gramigna 
GD, Lundy D, McGarvey-Toler S, Miller Gardner PJ: A randomized study of three interventions for aspiration of 
thin liquids in patients with dementia or Parkinson's disease. J Speech Lang Hear Res 2008, 51:1 73. 



Jestrovic et al. BioMedical Engineering OnLine 201 3, 1 2:90 
http://www.biomedical-engineering-online.eom/content/12/1/90 



Page 1 7 of 1 7 



23. Steele C, Van Lieshout P: Influence of bolus consistency on lingual behaviors in sequential swallowing. 

Dysphagia 2004, 1 9(3): 1 92-206. 

24. Pelletier C, Lawless H: Effect of citric acid and citric acid-sucrose mixtures on swallowing in neurogenic 
oropharyngeal dysphagia. Dysphagia 2003, 18(4):23 1-241. 

25. Lee J, Steele C, Chau T: Time and time-frequency characterization of dual-axis swallowing accelerometry 
signals. Physiol Meas 2008, 29(9):1 105-1 122. 

26. Amft O, Troster G: Methods for detection and classification of normal swallowing from muscle activation and 
sound. In Proc. of 1st International Conference on Pervasive Computing Technologies for Healthcare. Innsbruck, Austria; 
Nov.29-Dec. 1st, 2006:1-10. 

27. Reynolds EW, Vice FL, Gewolb IH: Variability of swallow-associated sounds in adults and infants. Dysphagia 

2009, 24:13-19. 

28. Shanahan TK, Logemann JA, Rademaker AW, Pauloski BR, Kahrilas PJ: Chin-down posture effect on aspiration in 
dysphagic patients. Arch Phys Med Rehabil 1 993, 74(7)736-739. 

29. Bulow M, Olsson R, Ekberg O: Supraglottic swallow, effortful swallow, and chin tuck did not alter 
hypopharyngeal intrabolus pressure in patients with pharyngeal dysfunction. Dysphagia 2002, 1 7(3):1 97-201 . 

30. Ertekin C, Keskin A, Kiylioglu N, Kirazli Y, On AY, Tarlaci S, Aydogdu I: The effect of head and neck positions on 
oropharyngeal swallowing: a clinical and electrophysiologic study. Arch Phys Med Rehabil 2001, 82(9)1 255. 

31. Bulow M, Olsson R, Ekberg O: Videomanometric analysis of supraglottic swallow, effortful swallow, and chin 
tuck in healthy volunteers. Dysphagia 1999, 14(2):67-72. 

32. Adnerhill I, Ekberg O, Groher ME: Determining normal bolus size for thin liquids. Dysphagia 1989,4:1-3. 

33. Sejdic E, Komisar V, Steele C, Chau T: Baseline characteristics of dual-axis cervical accelerometry signals. Ann 
BiomedEng 2010, 38(3): 1 048-1 059. 

34. Sejdic E, Steele C, Chau T: A method for removal of low frequency components associated with head 
movements from dual-axis swallowing accelerometry signals. PLoS One 201 2, 7(3):33464-33473. 

35. Donoho D: De-noising by soft-thresholding. IEEE Trans Inf Theory 1995, 41 (3):61 3-627. 

36. Sejdic E, Steele C, Chau T: Segmentation of dual-axis swallowing accelerometry signals in healthy subjects 
with analysis of anthropometric effects on duration of swallowing activities. IEEE Trans Biomed Eng 2009, 
56(4):1 090-1 097. 

37. Porta A, Guzzetti S, Montano N, Furlan R, Pagani M, Malliani A, Cerutti S: Entropy, entropy rate, and pattern 
classification as tools to typify complexity in short heart period variability series. IEEE Trans Biomed Eng 2001 , 
48(1 1):1 282-1 291. 

38. Porta A, Basel li G, Liberati D, Montano N, Cogliati C, Gnecchi-Ruscone T, Malliani A, Cerutti S: Measuring regularity 
by means of a corrected conditional entropy in sympathetic outflow. Biol Cybern 1 998, 78:71-78. 

39. Lem pel A, Ziv J: On the complexity of finite sequences. IEEE Trans Inf Theory 1 976, 22:75-81 . 

40. Rosso O, Blanco S, Yordanova J, Kolev V, Figliola A, Schurmann M, Basar E: Wavelet entropy: a new tool for 
analysis of short duration brain electrical signals. J Neurosci Methods 2001 , 1 05:65-76. 

41 . Pittner S, Kamarthi S: Feature extraction from wavelet coefficients for pattern recognition tasks. lEEETrans 
Pattern AnalMach Intel 1 1 999, 21 :83-88. 

42. Yen G, Lin K: Wavelet packet feature extraction for vibration monitoring. IEEE Trans Ind Electron 2000, 
47(3):650-667. 

43. DePuy V, Berger V, Zhou Y: Wilcoxon-Mann-Whitney test. In Encyclopedia of Statistics in Behavioral Science. Edited 
by Everitt B, Howell D. New York: John Wiley & Sons; 2005:21 18-2121. 

44. Lee J, Sejdic E, Steele C, Chau T: Effects of liquid stimuli on dual-axis swallowing accelerometry signals in a 
healthy population. Biomed Eng Online 201 0, 9:7-21 . 

45. Zhang X, Roy R, Jensen E: EEG complexity as a measure of depth of anesthesia for patients. lEEETrans Biomed 
Eng 2001, 48(1 2):1 424-1 433. 

46. Szczepanski J, Amigo J, Wajnryb E, Sanchez-Vives M: Application of Lempel-Ziv complexity to the analysis of 
neural discharges. Netw: Comput Neural Syst 2003, 14(2):335-350. 

47. Youmans S, Stierwalt J: An acoustic profile of normal swallowing. Dysphagia 2005, 20(3):1 95-209. 

48. Hamlet S, Patterson R, Fleming S, Jones L: Sounds of swallowing following total laryngectomy. Dysphagia 1 992, 
7(3):160-165. 

49. Reimers-Neils L, Logemann J, Larson C: Viscosity effects on EMG activity in normal swallow. Dysphagia 1 994, 
9(2):101-106. 

50. Dantas R, Kern M, Massey B, Dodds W, Kahrilas P, Brasseur J, Cook I, Lang I: Effect of swallowed bolus variables on 
oral and pharyngeal phases of swallowing. Am J Physiol-Gastrointest Liver Physiol 1 990, 258(5):G675-G681 . 

doi:1 0.1 1 86/1 475-925X-1 2-90 

Cite this article as: Jestrovic etal:. The effects of increased fluid viscosity on swallowing sounds in healthy adults. 

BioMedical Engineering OnLine 201 3 1 2:90. 
V J 



